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ABSTRACT 

The redshifted light of the first (Pop III) stars might substantially contribute to the 
near infrared background (NIRB). By fitting recent data with models including up- 
to-date Pop III stellar spectra, we find that such stars can indeed account for the 
whole NIRB residual (i.e. after 'normal' galaxy contribution subtraction) if the high 
redshift star formation efficiency is /* = 10% — 50%, depending on the IMF (the 
top-heaviest requiring lowest efficiency) and on the unknown galaxy contribution in 
the L band (our models, however, suggest it to be negligible). Such epoch of Pop 
III star formation ends in all models by z enc i w 8.8, with a hard limit z enc i < 9 set 
by J band observations. To prevent an associated IGM over-enrichment with heavy 
elements compared to observed levels in the IGM, pair-instability supernovae must be 
the dominant heavy element sources. Alternative explanations must break the light- 
metal production link by advocating very massive stars (M > 260M Q ) locking their 
nucleosynthetic products in the compact remnant or by postulating an extremely 
inhomogeneous metal enrichment of the Lya forest. We discuss these possibilities in 
detail along with the uncertainties related to the adopted zodiacal light model. 

Key words: galaxies: formation - intergalactic medium - black holes - cosmology: 
theory 



1 INTRODUCTION 

The cosmic infrared background (CIRB) can set important 
constraints on the star formation history of the universe (e.g. 
Hauser & Dwek 2001 and references therein). The CIRB 
consists of the cumulative radiation from all extragalactic 
sources. In the near-infrared (NIR) the CIRB is dominated 
by starlight, whereas in the mid- and far-infrared it results 
from dust emission. 

Different authors had provided measures of the Near 
Infrared Background (NIRB) in the J, K, and L bands us- 
ing the data of DIRBE instrument aboard of COBE satellite. 
Recently, Matsumoto et al. (2000) have accurately measured 
the NIRB from 1.4 to 4 /im, using the data taken with the 
Near Infrared Spectrometer. The largest uncertainty on the 
NIRB measurements is due to the subtraction of the in- 
terplanetary dust (IPD) scattered sunlight. In spite of this 
difficulty some studies have concluded that normal galaxies 
cannot account for the whole observed NIRB (Totani et al. 
2001). 

Bond, Carr & Hogan (1986) pointed out that the first 
population of zero-metallicity stars, the so called Popula- 
tion III (Pop III), might contribute to the cosmic back- 
ground in the near infrared. More recently, Santos, Bromm 



& Kamionkowski (2001) suggested that a significant part of 
the unaccounted NIRB could come from Pop III stars. They 
considered a very top-heavy Initial Mass Function (IMF), as- 
suming that all stars have masses > 300 Mq as recent three- 
dimensional numerical simulations seem to suggest (Bromm, 
Coppi & Larson 1999, 2002; Abel, Bryan & Norman 2000). 
In this paper we improve on these earlier studies in many 
respects. First, we confront theoretical predictions with the 
recent Matsumoto et al. (2000) data, which provides a much 
tighter test and constraint. Secondly, we consider the effect 
of several IMFs. Third, we use the new stellar spectra of 
Schaerer (2002) for zero age mean sequence (ZAMS) metal 
free stars, including the nebular continuous emission found 
to be very important for stars with strong ionizing fluxes. 

We fit the estimate of the unaccounted NIRB finding 
limits on the star formation efficiency and on the redshift 
at which the formation of Pop III stars ends. We compare 
our results with high redshift observations and draw some 
results on IGM metal enrichment due to Pop III supernovae. 

The paper is organized as follows: in §2 we describe the 
background intensity calculation, the stellar spectra used, 
and how we model the IGM. In §3 the available data in 
the NIR range are reviewed along with the number counts 
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of normal galaxies from deep field surveys. We discuss the 
results of our analysis in §4; in §5 we present our conclusions. 

We adopt the 'concordance' model values for the cos- 
mological parameters: h — 0.7, fi M = 0.3, Qa = 0.7, 
n b = 0.038, ct 8 = 0.9, and T = 0.21. Here h is the di- 
mensionless Hubble constant, Ho = 100/t km s _1 Mpc -1 ; 
Q.M, Oa, and Q b are the total matter, cosmological con- 
stant, and baryon density in units of the critical density; erg 
gives the normalization of the power spectrum on 8/i -1 Mpc 
scale and T is the shape of the power spectrum. We use the 
power spectrum for the fluctuations derived by Efstathiou, 
Bond & White (1992). 



2 COSMIC INFRARED BACKGROUND 

The mean specific intensity of the background 7(Vo,zo), as 
seen at frequency by an observer at redshift zo, is given 



by 



/(«*,*>) = ^ [°° e{v,z)e- T >ff<™-*^dz 

J zn 



(1) 



(Peebles 1993). Here e(y,z) is the comoving specific emis- 
sivity, v — fo(l + z)/(l + zo), T e ff(vo, zo, z) is the effective 
optical depth at of the IGM between redshift zo and z, 
and 61/ Az is the proper line element 



^ z =c[H (l + z)E(z)]- 1 , 
where c is the light speed and 

e(z) = [n M (i + zf + n A + (i - n M - n A )(i + zf] 1 ' 2 . (2) 



2.1 The Comoving Specific Emissivity 

The gas initially virialized in the potential well of the par- 
ent dark matter halo, can subsequently fragment and ignite 
star formation only if the gas can cool efficiently and lose 
pressure support. For a plasma of primordial composition 
at temperature T < 10 4 K, the typical virial temperatures 
of the early bound structures, the only efficient coolant is 
molecular hydrogen. Thus, a minimum H2 fraction is re- 
quired for a gas cloud to be able to cool in a Hubble time. 
As the intergalactic relic H2 abundance falls short by at 
least two orders of magnitude such requirement, the fate of 
a virialized lump depends crucially on its ability to rapidly 
increase its H2 content during the collapse phase. Once a 
critical H2 fractional abundance of ~ 5 x 10~ 4 is achieved 
in an object, the lump will cool, fragment and eventually 
form stars. This criterion is met only by larger haloes, so 
that for each virialization redshift there will exist a criti- 
cal mass, M m i n (z), for which protogalaxies with total mass 
Mh > M m in form stars and those with mass Mh < M m i n 
do not. We adopt here M m i„(z) as computed by Fuller & 
Couchman (2000). 

In absence of additional effects that could prevent or 
delay the collapse, we can associate to each dark matter 
halo with mass M h > M min a corresponding stellar mass 



\Im 



where /* is the fraction of baryons able to cool and form 
stars; we will refer to this quantity as the star formation 
efficiency. Thus, the stellar mass per unit comoving volume 
at redshift z contained in haloes with mass M h > M min (z) 
is given by 



POO 

p*(z)= / n(M h ,z)M*AM h , 

J M min {z) 



(4) 



where n(M h ,z) is the comoving number density of dark 
matter haloes of mass M h at redshift z given by Press & 
Schechter (1974). 

The comoving specific emissivity, in units of erg s _1 
Hz -1 cm~ 3 , is then 



e(v,z) = l v (z)p+(z), 



(5) 



where l„(z) is the specific luminosity of the population [in 
erg s -1 Hz -1 Mg 1 ] at redshift z (see Sec. 2.4). 



2.2 The Intergalactic Medium 

Absorption due to intergalactic gas located in discrete sys- 
tems along the line of sight can seriously distort our view 
of objects at cosmological distances. At wavelengths short- 
ward of Lya (A < 1216 A) in the emitter rest frame, the 
source continuum intensity is attenuated by the combined 
blanketing of lines in the Lyman series, and strongly sup- 
pressed by the continuum absorption from neutral hydrogen 
shortward of the Lyman limit in the emitter rest frame. As 
pointed out by Haiman & Loeb (1999) no flux is transmit- 
ted at wavelengths Ao < A a (l + z s ) from sources at redshift 
1 + z s > 32/27(1 + Zi), where z; = 6.2(Gnedin 2001) is the 
currently favored reionization redshift. 

The effective optical depth r e // through the IGM is 
defined as e~ Tcff —< e~ T >, where the mean is taken over 
all the lines of sight to the redshift of interest. For a Poisson 
distribution of absorbers (Madau 1991, 1992), 



f-Z POO 

T eff {vo,zo,z) = dz dN Rl C(N Rl ,z')(l-e~ T ), 

J zn JO 



(6) 



where C(-^Hi > z ') ~ d 2 N/dN-Q 1 dz' is the distribution of 
the absorbers as a function of redshift and neutral hydro- 
gen column density, JVjj r ; r(z/) is the optical depth of an 
individual cloud for ionizing radiation at frequency v 1 . 

For the Lyman-a forest, in the wavelength range A/3 < 
A /(l + z em ) < \ a where X a = 1216 A and A^ = 1026 A, 
we have (Madau 1995) 



0.0036 



(aJ 



When Ao/(l + z em ) < A/3, a significant contribution to 
the blanketing opacity comes from the higher order lines of 
the Lyman series. In the wavelength range A»+i < Ao/(l + 
z) < Aj, the total effective line-blanketing optical depth can 
be written as the sum of the contributions from the j — > 1 
transitions, 



(3) 



Teff 



(7) 
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= 1.5 x 10-Hz ( <h*) (—) " ii±^!, (11) 
n ' ^0.019 J \0.7J E{z) ' v ; 

where E(z) is defined in eq. (2). Note that we have allowed 
the IGM in the vicinity of the emitting galaxy to be over- 
dense by a factor 6 « 10, to account for the high clustering 
of first sources. 



Table 1. Coefficients Aj corresponding to Xj in eq. 7 



2.4 Emission from a Pop III 'Stellar Cluster' 
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Table 2. Best fit value (Model Al) for A, f3 and 7 (Fardal et al. 
1998) used in eq. 9 



where Aj and the corresponding values for Aj are given in 
Table 1. 

Continuum absorption from neutral hydrogen along the 
line of sight affects photons observed at Ao/(l + z em ) < Al, 
where Al = 912 A is the Lyman limit. We have 



Hi a 

where a(X ,z) ~ 6.3 x KT 18 (A /A L ) 3 (1 + z)~ 



(8) 

ti 2 (Oster- 

brock 1989) is the hydrogen photoionization cross section, 
and (1 + z c ) = (Ao/Al) for Ao > Al, whereas for Ao < Al, 
Hi absorbs photons all the way down to z c — 0. 

For the redshift and column density distribution of ab- 
sorption lines, the usual form can be adopted 



C(^Hi ,*) 



Go 17 ) 



Hi 



10 17 cm - 



{1 + zf 



(9) 



where the values of coefficients A, (3 and 7 in different ranges 
in iVjjj are given in Table 2, and are taken from Fardal, 
Giroux & Shull (1998). 



2.3 Scattering of Lya Photons 

Lya line photons from first galaxies are absorbed by neutral 
hydrogen along the line of sight (Gunn & Peterson 1965). 
However, they are not destroyed but scatter and diffuse in 
frequency to the red of the Lya resonance owing to the 
Hubble expansion of the surrounding Hi . Eventually, when 
their net frequency shift is sufficiently large, they escape 
and travel freely toward the observer. The profile results in 
a strong, asymmetric Lya emission line around 1225 A with 
a scattering tail extending to long wavelengths. 

The resulting scattered line profile, $(v, z), has been 
simulated by Loeb & Rybicki (1999). Here, we use the fitting 
analytical expression given by Santos et al. (2001) 



u i ,(z)iy- 2 exp ifi/>0 
if v < 



(10) 



We calculate the emission of a Pop HI 'stellar cluster' ac- 
cording to the equation: 



l v (z) = / F(v,M,z)<t>(M)AM, 

J Mt 



(12) 



where 4>(M) is the IMF normalized so that f™ u (p(M)dM = 
1; Mi and M u are the lower and upper mass limits. The 
spectrum of a star of mass M at the rest-frame frequency v 
and redshift z is given by 



F(v, M, z) = C ar (M) + IT°(M) + t ya (M, z). 



(13) 



where ll tar (M) is the spectrum of a Pop III star of mass 
M, l™ eb is the emission of the nebula surrounding the star, 
and li yoi {z) is the emission due to Lya line photons from a 
source at redshift z scattered by the IGM (see Sec. 2.3). 

As spectra of Pop HI we adopt here those modeled by 
Schaerer (2002) who relaxed some assumptions made in pre- 
vious works (Tumlinson & Shull 2000; Tumlinson, Giroux 
& Shull 2001; Bromm, Kudritzki, & Loeb 2001). Schaerer 
(2002) presented realistic models for massive Pop III stars 
and stellar populations based on non-LTE atmospheres, re- 
cent stellar evolution tracks and up-to-date evolutionary 
synthesis models. Moreover, he included nebular continuous 
emission, which cannot be neglected for metal-poor objects 
with strong ionizing fluxes. This process increases signifi- 
cantly the total continuum flux at wavelengths redward of 
Lyman-a and leads in turn to reduced emission line equiv- 
alent widths. 

We add thus the contribution of the nebular emission 
(n e = 100 cm -3 is assumed) as described in Schaerer (2002): 



r 
'-I' 



OLB 



(14) 



(e.g. Osterbrock 1989), where as (in units of cm~ 3 s" 1 ) 
is the Case B recombination coefficient for hydrogen and 
/ esc is the ionizing photon escape fraction out of the ideal- 
ized region considered here (we assume f esc — 0); q(H) = 
Q(H)/M, where Q(H) is the ionizing photon rate (in units 
of s _1 ) for H and is given for different stellar masses in Ta- 
ble 3 of Schaerer (2002). The continuous emission coefficient 
ftot, including free-free and free-bound emission by H, neu- 
tral He and single ionized He, as well as the two-photon 
continuum of hydrogen is given by 



7tot = IHI + "/2q + iHel 



n{He 



+ Well- 



(15) 



The continuous emission coefficients 7; (in units of erg 
cm 3 s _1 Hz -1 ) are taken from Aller (1987) for an electron 
temperature of 20000 K. The spectrum of a 1000 Mq star 
(with and without the nebular emission) and of a 5 M Q star 
are plotted in Fig. 1. 
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hu\ eV] 

Figure 1. Spectra of individual Pop III stars plotted as luminos- 
ity per unit of stellar mass (in erg s _1 Hz -1 M^ 1 ) vs. energy (in 
eV). The spectrum of a 1000 Mq star with (solid line) and with- 
out (dashed) nebular emission and of a 5 Mq are shown. For the 
5 Mq star the contribution of the nebular emission is negligible. 



Finally, the emission of Lya line photons is given by 
t ya (z) = c Lya (l - f esc ) q(H) <b(v Lya -v,z), (16) 



where c Ly c 



1.04 x 10" 



erg (Schaerer 2002). The line 



profile $(^) is given in eq. (10). 
2.4-1 Initial Mass Function 

The issue of the mass distribution of first metal free stars 
has been tackled via various hydrodynamical models and 
other studies (e.g. Abel et al. 1998, Bromm et al. 1999, 2001, 
2002, Nakamura & Umemura 2001, Ripamonti et al. 2002). 
There seems to be an overall consensus that massive stars 
(up to 1000 Mq) may form. Simulations seem to suggest 
that the primordial IMF might have been biased towards 
stellar masses > 100 Mq, but other studies (e.g. Nakamura 
& Umemura 2001) found that the formation of stars with 
masses down to 1 Mq is not excluded. 

In the local universe the IMF is usually well described 
by the standard Salpeter law (Salpeter 1955) 

-2.35 



4>(M) oc M~ 



(17) 



Larson (1998) suggested a different form for the high z 
IMF. This has a universal Salpeter-like form at the upper 
end, but flattens below a characteristic mass, M c which may 
vary with time. By increasing M c we can mimick a top-heavy 
IMF of Pop III stars, 



^.(M)cxM- 1 (l + H) 



(18) 



Hernandez & Ferrara (2001) have explored the pre- 
dictions of the standard hierarchical clustering scenario of 
galaxy formation, regarding the numbers and metallicities 
of Pop III stars that are likely to be found within our Galaxy 
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1 20 



= 10 




10 
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Figure 2. Pop III stellar cluster spectra for different IMFs. 
Dashed line: Salpeter. Dotted line: heavy IMF. Solid line: very 
heavy IMF. The sharp peak around 1220 A is due to Lya nebu- 
lar emission for sources at z s = 10, filtered through the IGM. The 
flux at E> 13.6 eV is then strongly suppressed by the absorption 
due to the intcrgalactic gas (see Sec. 2.2) 



today. By comparing these values with observational data, 
they suggested that the IMF of first stars was increasingly 
high-mass weighted towards high redshifts, levelling off at 
z > 9 at a characteristic stellar mass scale of 10-15 Mq. 

In view of our poor knowledge of the IMF at high red- 
shift we consider here three different distributions of mass for 
the first population of stars: Salpeter, Larson with M c = 15 
Mq (heavy IMF), and Larson with M c = 100 Mq (very 
heavy IMF). For all IMFs, Mi = 1 Mq and M u = 1000 Mq. 
In Fig. 2 are shown the resulting spectra of a Pop III cluster 
at redshift 2 S =10 for the three different IMFs assumed. 



3 OBSERVATIONAL CONSTRAINTS 

Observations of the NIRB are seriously hampered by the 
strong atmospheric foreground (e.g. Mandolesi et al. 1998). 
On the other hand space measurements are not easy be- 
cause it is difficult to subtract correctly the contribution to 
the extragalactic background light by interplanetary dust 
scattered sunlight (zodiacal light). Using the data from the 
instrument DIRBE aboard of the COBE satellite and from 
the 2MASS data, Cambresy et al. (2001) found a non zero, 
isotropic background in the J and K bands. With the same 
data but a different zodiacal light model, Wright & Johnson 
(2001) extended the flux estimate to the L band. Cambresy's 
and Wright's results in the J and K bands are compatible if 
the same zodiacal light model is applied. Matsumoto et al. 
(2000) provided a preliminary analysis of the the NIRS data 
obtaining a detection of the NIRB in the 1.4-4 /im range. 
They applied the same zodiacal light model of Cambresy et 
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Figure 3. Available NIR data. The small filled circles are the 
NIRS data (Matsumoto et al. 2000). The open symbols are the 
DIRBE results: squares for Wright (2001), diamonds for Cam- 
bresy et al. (2001), and circles for Gorijan, Wright & Chary 
(2000). The big filled circle is the Kiso star count measurement. 
The data arc slightly offset for clarity. The errors are at la and for 
all the data the Kelsall et al. (1998) model for the zodiacal light 
is applied. The open triangles are the count integration from the 
Hubble Deep Field (Madau & Pozzetti 2000), whereas the filled 
triangles report upper and lower limits on the count integration 
from Subaru Deep Field when also the contribution of missed 
galaxies is considered (Totani et al. 2001). 



al. (2001) and obtained compatible results in the DIRBE 
bands. At 8000 A, Bernstein, Freedman & Madore (2002) 
observed a CIRB of 1.76 ± 0.48 x 10~ 5 erg s _1 cm~ 2 sr _1 , 
that is ~4 times smaller than the DIRBE and NIRS mea- 
surements in the J band, showing a break around ~1 /im. 
The integrated counts in the J and K bands from deep sur- 
veys do not account for the total NIRB (Madau & Pozzetti 
2000; Totani et al. 2001). Totani et al. (2001) modeled the 
contribution of galaxies missed by deep galaxy surveys. They 
found unlikely that the contribution of all the normal galax- 
ies to the NIRB is larger than 30%. All the available data 
are plotted in Fig. 3. In Table 3 the estimate of the NIRB 
from the DIRBE data (top panel) and the contribution of 
normal galaxies (bottom panel) are reported. 



3.1 The DIRBE Data 

The Diffuse Infrared Background Experiment (DIRBE) on 
board of the Cosmic Background Explorer (COBE, see 
Boggess et al. (1992)) satellite was designed to search for 
the Cosmic Infrared Background from 1.25 /im to 240 /im. 
The DIRBE instrument (Silverberg et al. 1993) was an ab- 
solute photometer which provided maps of the full sky in 
10 broad bands at 1.25, 2.2, 3.5, 4.9, 12, 25, 60, 100, 140, 
and 240 fim with a starlight rejection of < 1 nW m~ 2 sr _1 
and an absolute brightness calibration uncertainty of 0.05 
and 0.03 nW m~ 2 sr _1 , at 1.25 /im and 2.2 /j,m respectively. 



A summary of the DIRBE results is provided by Hauser et 
al. (1998) and all the NIRB measures based on the DIRBE 
data are reported in Table 3 (top panel). 

In the near infrared bands the dominant foreground in- 
tensities in the DIRBE data are the zodiacal light and the 
light from stars in the Milky Way. 

A first attempt to detect the NIRB from the DIRBE 
data was carried out by the DIRBE group but this gave 
only an upper limit and failed the tests of isotropy even in 
limited regions of sky (Arendt et al. 1998). Dwek & Arendt 
(1998) made a correlation study for the K and L bands and 
obtained a lower limit for the L band. 

Gorjian et al. (2000) removed the foreground due to 
Galactic stars by directly measuring all stars brighter than 
9th magnitude at 2.2 and 3.5 /jm in a 2° x 2° dark spot 
near the North Galactic Pole using ground-based telescopes. 
They calculated the contribution of fainter stars using the 
statistical model of Wainscoat et al. (1992) and subtracted 
the zodiacal light contribution using a improvement of the 
Wright (1998) model. Gorjian et al. found significant pos- 
itive residuals in the K and L bands which they identi- 
fied as a probable detection of the NIRB. Wright & Reese 

(2000) obtained a consistent estimate of the NIRB in the 
same bands using a different approach based on a histogram 
fitting method to remove the stellar foreground from the 
DIRBE data. 

Wright (2001) used the Two Micron All Sky Survey 
(2MASS) data (Cutri et al. 2000) to remove the contribution 
of Galactic stars brighter than 14th magnitude from the 
DIRBE maps at 1.25 and 2.2 fj,m in four dark regions in the 
north and south Galactic Pole caps. For the subtraction of 
the zodiacal light foreground the model presented in Gorjian 
et al. (2000) is applied. Cambresy et al. (2000) have also 
used the 2MASS data to remove the Galactic stars from 
the DIRBE data, but they modeled the zodiacal light as in 
Kelsall et al. (1998). Their results agree with those of Wright 

(2001) if the same zodiacal light model is applied. 
Recently, Wright & Johnson (2001) extended the analy- 
sis of Wright (2001) to 13 fields with a wide range of ecliptic 
latitudes tripling the number of pixels. They obtained an 
estimate of the NIRB at 1.25 and 2.2 fim consistent with 
the previous results. They also found a significant resid- 
ual in the L band combining the Wright (2001) and Dwek 
& Arendt (1998) techniques, consistent with the Wright & 
Reese (2000) result. 

Matsumoto et al. (in preparation) observed in the J 
band with the Kiso Schmidt Telescope toward the DIRBE 
dark spot. They subtracted all detected stars brighter than 
14th magnitude and estimated the contribution of fainter 
stars by the Cohen's sky model (Cohen 1997). They used 
the KelsalPs zodiacal light model. 



3.2 The NIRS Data 

The Near Infrared Spectrometer (NIRS)* is one of the fo- 
cal instruments of the InfraRed Telescope in Space (IRTS) 
(Noda et al. 1994). The NIRS covers the wavelength range 
from 1.4 to 4.0 fim with a spectral resolution of 0.13 /jm. 
The beam size is 8' x 20'. which is considerably smaller than 



http://www.ir.isas.ac.jp/irts/nirs/index-e.html 
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Wright & Reese 2000 


(27.7±14.5) 


(19.9±5.3) 
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Wright 2001 


60.8±14.5 


28.2±5.5 




K 




54.0±16.8 


27.8±6.7 
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Cambrcsy ct al. 2001 


(24.3±18) 


(24±6) 


(13.8±3.4) 


W 


Wright & Johnson 2001 



60.1±15 K Kiso star counts 



Simple Integration of Observed Galaxy Counts 



(9-7t?'i)° 7.9tl X 2 Madau & Pozzctti 2000 

10.9 ± 1.1 8.3 ± 0.8 Totani et al. 2001 



Estimated Resolved Fraction 



0.97 0.93 Model A 

0.95 0.92 Model B 



Contribution to the NIRB of All Normal Galaxy 
10.1-12.8 7.8-10.2 Totani et al. 2001 



Table 3. Summary of the recent observations of the NIRB from the DIRBE data (top panel). Units are nWm -2 sr _1 (or 10 -6 erg s _1 
cm 2 sr -1 ). In the fourth column is given the Model used to subtract the Zodiacal Light: 'W stands for Wright (1998) and 'K' is for 
Kclsall ct al. (1998). In the bottom panel are shown the results of count integration from the Hubble Deep Field (Madau & Pozzctti 
2000) and from the Subaru Deep Field (Totani et al. 2001) surveys, the estimate of the contribution of missed galaxies (Totani et al. 
2001. Model A: number evolution with n ~ 1; Model B: no number evolution with r\ = 0) and the total flux from all normal galaxies 
(Totani ct al. 2001). 

a Estimate at slightly different wavelengths 



that of DIRBE; about 7% of the sky is surveyed. In order to 
reduce the contribution from the faint stars, the sky at high 
latitudes (6 > 40°) is chosen. 

Matsumoto et al. (2000) provided a preliminary anal- 
ysis of the NIRS data, reporting detection of the Cosmic 
Infrared Background based upon analysis of the five days of 
data which were least disturbed by atmospheric, lunar, and 
nuclear radiation effects. The sky area analyzed included 
Galactic latitudes from 40° to 58°, and ecliptic latitudes 
from 12° to 71°. NIRS is able to identify stars brighter than 
10.5 mag. at 2.24 /u,m, whereas to find out the contribution 
of fainter stars the Cohen model (1997) is used. The model 
of Kelsall et al. (1998) is applied to subtract the contribu- 
tion of the zodiacal light foreground, interpolating between 
the DIRBE wavelengths to the wavelengths of the NIRS 
measurements. After subtraction of the IPD contribution, 
there remains a fairly isotropic residual emission, which they 
interpreted as evidence for a non-zero background. To ob- 
tain a quantitative value for the background at each wave- 
length they correlated their star-subtracted brightness at 
each point with the IPD model brightness, and used the ex- 
trapolation to zero IPD contribution as a measurement of 
the NIRB. 

The NIRB intensities reported by Matsumoto et al. 
(2000) near 2.2 and 3.5 jj,m are similar to the values found 
by Gorjian et al. (2000). At shorter wavelengths, the NIRS 



results continue to rise steeply to ~ 6.5 x 10~ 5 erg s _1 cm~ 2 
sr -1 at 1.4 /im. This is somewhat above the 95% CL upper 
limit at 1.25 fim of Wright (2001), but it is in agreement 
with the value obtained by Cambresy et al. (2001) using 
the same zodiacal light model. The NIRS result implies an 
integrated background energy over the 1.4-4.0 /itm range of 
~ 3.0 x 10~ 5 erg s _1 cm~ 2 sr _1 . 

The preliminary report by Matsumoto et al. (2000) does 
not provide details regarding systematic uncertainties in 
their results. The largest uncertainty in the reported NIRB 
values is attributed to the IPD model. 



3.3 Galaxy Contribution from Deep Counts 

Deep optical and NIR galaxy counts provide an estimate 
of the extragalactic background light (EBL) coming from 
normal galaxies in the universe. Madau & Pozzetti (2000) 
derived the contribution of known galaxies in the UBVI- 
JHK bands from the Southern Hubble Deep Field. Although 
the slope of the number-magnitude relation of the faintest 
counts is flat enough for the count integration to converge, 
the recent observations of the NIRB (see previous Section) 
suggest that the diffuse EBL flux is considerably larger than 
the count integration. However, a considerable fraction of 
EBL from galaxies could still have been missed in deep 
galaxy surveys because of various selection effects. 
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flat 


0.24±0.03 


0.39±0.05 


0.58±0.08 


Heavy 
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0.15±0.01 


0.24±0.02 


0.36±0.03 
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0.15±0.02 
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0.10±0.01 
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0.15±0.02 



Table 4. Values of /* from the fit to the data with different IMFs. 
The quoted statistical errors are at 95% C.L.. 
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Totani et al. (2001) using the Subaru Deep Field data 
(Maihara et al. 2001) have shown that more than 80%-90% 
of EBL from galaxies has been resolved in the K and J bands 
and that the contribution by missing galaxies cannot ac- 
count for the discrepancy between the count integration and 
the NIRB observations. In Table 3 (bottom panel) we show 
the count integration in optical and NIR bands from Madau 
& Pozzetti (2000) and Totani et al. (2001), and the resolved 
fractions for two galaxy evolution models. The possible num- 
ber evolution of galaxies is considered by a phenomenological 
model in which the Schechter parameters of the luminosity 
function have a redshift dependence as (f>* oc (1 + z) v and 
L* oc (1 + z)~ v ; i.e., luminosity density is conserved. In 
Model A an evolution with ij = 1 is considered, whereas in 
Model B no evolution (rj = 0) is assumed (see Totani et al. 
(2001) for more details). Currently no measure of the contri- 
bution of normal galaxies to the NIRB in L band is available 
whose accurate measurement has to await for SIRTF. 



4 MODEL VS. OBSERVATIONS 

We compare the NIR data with the predictions of our model, 
summarized by eq. (1). The data set is composed by the mea- 
surements of Matsumoto et al. (2000) and those obtained 
by the DIRBE data using the Kelsall model for the zodiacal 
light emission as reported in Table 3 (top panel). For each 
wavelength we subtracted from the data the contribution 
of normal galaxies using the upper (lower) limits in the J 
and K band, obtained by Totani et al. (2001) (last line of 
Table 3). As discussed already, for A > 2.2 /jm no estimate 
of the NIRB due to normal galaxies is available. Thus in 
the range 2.2 to 4 fim we consider two extreme cases: no 
contribution from normal galaxies (labeled with cut) and a 
constant contribution (labeled with flat). 

4.1 Model Results 

The fit results are quite sensitive to the adopted IMF and 
to the way we model the contribution of normal galaxies for 
A > 2.2 fim. In Table 4 are reported the results of the fit 
for different IMFs, count integration model, and z en d- The 
reported errors represent the statistical uncertainty of the 
fit, but systematic uncertainties difficult to evaluate might 
be present. In Fig. 4 are plotted the values of the star for- 
mation efficiency and the corresponding x 2 /DOF (Degree of 
Freedom) as function z en d for the cut galaxy count model. 




8 8.5 9 9.5 10 



Z end 

Figure 4. Star formation efficiency /* accounting for the whole 
unexplained NIRB, as a function of redshift at which the forma- 
tion of Pop III stars ends. The model cut for the counts inte- 
gration for A > 2.2 (im is adopted. Dashed line: Salpeter IMF. 
Dotted line: heavy IMF. Solid line: very heavy IMF. The bottom 
panel shows the value of x 2 /DOF. 



The x 2 increases rapidly as z en d becomes greater than 9, 
showing that in this case our model fails to account of the 
NIRB observation in the J band. 

For a Salpeter IMF we find that a large fraction of 
baryons have to be turned into stars to account of the whole 
unexplained NIRB. This fraction goes from 37% (z en d = 8) 
to 89% (z en£ j = 10) if we consider the case of no contribution 
from normal galaxies for A > 2.2 /im, slightly dependent on 
the choice of the upper and lower limits of the count inte- 
gration. If we consider a constant (flat model) contribution 
of normal galaxies up to 2.2 /mi the value of /* goes from 
~ 0.24 (z end = 8) to ~ 0.58 (z end = 10). The value of /* de- 
creases considerably if we consider a top-heavy IMF. For the 
cut model of the normal galaxy contribution, /* ranges be- 
tween ~ 0.15 (z en d = 8) and ~ 0.36 (z en d = 10) if M c — 15 
M Q and between ~ 0.1 and ~ 0.23 if M c = 100 Mq , whereas 
for the flat model we find /* values from ~ 0.1 (z end = 8) 
to ~ 0.23 (z end = 10) with a heavy IMF, and ~ 0.06 - 0.15 
(for z en d = 8 — 10) with a very heavy IMF. 

To find the best model, we allow both /* and z end to 
vary simultaneously and perform a two-parameter fit to the 
NIR data. The results are reported in Table 5. Again, quoted 
errors take into account for statistical uncertainties of the 
fit. In Fig. 5 we plot the best fit models for the cut contribu- 
tion of normal galaxies. For the flat model the star forma- 
tion efficiencies are significantly lower than in the cut case, 
but the x 2 /DOF remains always greater than 1.9, showing 
that a low contribution of normal galaxies for A > 2.2 /im is 
favored by our model. A very remarkable result is the con- 
stancy of z end fa 8.8 in the various cases, independent on the 
IMF or the galaxy contribution model: its value is strongly 
constrained by the position of the break around 1 /jm. The 
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cut 
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8.79±0.10 




flat 


0.32±0.04 


8.69±0.01 
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cut 


0.21±0.02 


8.80±0.01 




flat 


0.13±0.02 


8.75±0.01 


Very 


cut 


0.14±0.01 


8.83±0.01 


heavy 


flat 


0.09±0.01 


8.79±0.01 



Table 5. Best fit values for /* and z m d with different IMFs. The 
quoted statistical errors are at 95% C.L.. 



D.0001 




Figure 5. Best results of the fit for different IMFs. Dashed line: 
Salpeter IMF z m d = 8.79, /* = 0.53. Dotted line: heavy IMF 
z end = 8.8, fi, = 0.21 . Solid line: very heavy IMF z end = 8.82, 
/* = 0.14. The model cut for the counts integration for A > 2.2 
(im is adopted. The data are explained in Sec. 3 and in Fig. 3 



break is principally due to the Lya line emission redshifted 
into the J band, so that A (, s ~ ^Ly a /(1 + Zend)- 

Barkana (2002) derived a value for /* using the mea- 
sured distribution of star formation rates in galaxies at vari- 
ous redshifts (Lanzetta et al. 2002) and semi-analytic models 
of hierarchical galaxy formation in a ACDM cosmology. He 
found a best fit value of 2.3lg'g% (2a), neglecting dust ob- 
scuration, but the uncertainties on this evaluation are still 
large. If this value applies also to the first star formation, 
for a Salpeter IMF the Pop III contribution to the NIRB is 
not significant (< 10%), whereas Pop III stars characterized 
by a top-heavy IMF continue to contribute substantially to 
the NIRB. 



4.2 Associated Metal Enrichment 

Recent theoretical analyses on the evolution of metal-free 
stars predict that the fate of the massive metal-free stars 



can be classified as follows (Heger, Woosley & Waters 2000, 
Chiosi 2000, Heger & Woosley 2002): 

(i) M > 260 M : the nuclear energy release from the 
collapse of stars in this mass range is insufficient to reverse 
the implosion. The final result is a very massive black hole 
(VMBH) locking up all heavy elements produced; 

(ii) 130 M Q < M < 260 M Q : the mass regime of the 
pair- instability supernovae (SN 77 ). Precollapse winds and 
pulsations should result in little mass loss, the star implodes 
reaching a maximum temperature that depends on its mass 
and then explodes, leaving no remnant. The explosion expels 
metals into the surrounding ambient ISM; 

(iii) 40 M Q < M < 130 M : black hole formation is the 
most likely outcome, because either a successful outgoing 
shock fails to occur or the shock is so weak that the fall- 
back converts the neutron star remnant into a black hole 
(Fryer 1999); 

(iv) 8 M Q < M < 40 M Q : SN explosion from 'normal' 
progenitors. 

Stars in the mass range (i) and (iii) above fail to eject 
most of (or all) their heavy elements, whereas (ii) and (iv) 
type of stars eject metals into the IGM. 

To estimate the metal production of Pop III objects 
we consider SN yields of Woosley & Weaver (1995) for pro- 
genitor metal-free stars between 8 and 40 M and the re- 
sults of Heger & Woosley (2002) for very massive objects 
of pair instability SNe originating from stars in the mass 
range ~ 130M Q to 260 Mq. We neglect the contribution 
from long lived intermediate mass stars (1 < M < 8 Mq) 
as their evolution time scale is longer than the Hubble time 
at the relevant redshift. The density of metals ejected into 
the IGM in units of the critical density, p c = 3Hq/8ttG, is 
obtained by 



fiz(z) = fz(<t>) 



(19) 



where fz is the fraction of metal ejected from normal and 
pair-instability supernovae; it depends on the adopted IMF, 
<f>. Finally, fi*(z) = p*/p c , where p* is given in eq. (4). 

We calculate the IGM metallicity by imposing that the 
whole unaccounted NIRB is due to Pop III stars. We use 
the best results for the efficiency reported in Fig. 5 for the 
cut model of the galaxy counts. We assume here that all 
the metal escape into the IGM and that they are smoothly 
distributed. Fig. 6 shows the metallicity of the IGM (in solar 
units) as a function of redshift. 

We compare our results with the estimated value of met- 
als in the Lya forest (10 14 5 < JV Rl < 10 16 ' 5 cm" 2 ) at z ~ 3 
(Songaila & Cowie 1996; Dave et al. 1998) and in the 'true' 
IGM (N m < 10 14 cm" 2 ) at z ~ 5 (Songaila 2002). Iden- 
tification of C IV, Si IV and O VI absorption lines which 
correspond to Lya absorption lines in the spectra of high 
redshift quasars has revealed that the low density IGM has 
been enriched up to Zigm ~ 

10 -2.5±o.5 Zq _ For lower col _ 

umn density clouds (the 'true' IGM) Songaila (2002) found 
a lower limit on the IGM metallicity at z ~ 5 of 10~ 3 ' 5 Zq. 

With the assumed star formation efficiency (/* = 0.534 
for a Salpeter; /* = 0.214 for a heavy IMF; /* = 0.142 for a 
very heavy IMF) the IGM would be enriched to the observed 
metallicity already at z = 15 — 18 for the Lya forest value 
and at z = 20 - 25 for the 'true' IGM lower limit (Fig. 6). 
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Figure 6. IGM mctallicity as a function of redshift. Dashed line: 
Salpcter IMF. Dotted line: heavy IMF. Solid line: very heavy 
IMF. The curves are obtained imposing that the whole unac- 
counted NIRB is due to Pop III stars and the best results for star 
formation efficiency of Fig. 5 are assumed. The horizontal solid 
lines represent the value of the metallicity in the Lya forest (Dave 
et al. 1998) and in the 'true' IGM (Songaila 2002). 



As seen from Fig. 6, for z ~ 8.8 the mean IGM metallicity 
exceeds by one order of magnitude that observed in the Lya 
forest. 

As an alternative, we are thus forced to consider the 
analogous case in which only pair-instability SNe eject met- 
als into the IGM (Fig. 7) . For a Salpeter IMF the obtained 
metallicity is now consistent with the limits of the Lya for- 
est, whereas for a heavy IMF this value is slightly exceeded. 
For top heavy IMFs, for which the contribution of normal 
SNe is not so important, the results change only slightly. 

In order to further reduce these metallicity values, we 
have to consider lower star formation efficiencies. Conse- 
quently only a fraction of the NIRB can be due to Pop III 
stars. By imposing that Pop III stars do not eject more met- 
als than the observed ones, we find self consistent limits on 
the star formation efficiency and on the contribution of the 
first population of stars to the NIRB; Table 6 summarizes 
the results. Apart for the case of a Salpeter IMF in which 
only pair-instability SNe contribute to the metal enrichment 
of the IGM, the star formation efficiency has to be small 
(< 10%) and the fraction of NIRB due to Pop III stars is 
< 20%. 

The value obtained by Songaila (2002) for the 'true' 
IGM at z ~ 5 allows us to set a lower limit on the star 
formation efficiency and the contribution of Pop III stars to 
the NIRB. The efficiency has to be greater than fewxl0~ 4 
to enrich the IGM to Z ~ 10~ 3 ' 5 Z Q ; thus, the Pop III 
star contribution to the NIRB is at least of the order of 
0.4-1.7%. This fraction increases up to 8% if we consider 
the case in which only pair-instability SNe contribute to the 
metal enrichment of the universe. 
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Figure 7. As Fig. 6 but in the case in which only pair-instability 
SNe eject metals into the IGM. Dashed line: Salpeter IMF. Dotted 
line: heavy IMF. Solid line: very heavy IMF. 
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3.8 x 10~ 2 
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0.6 x 10~ 2 


4% 




only SNj~/ 
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6% 



Table 6. The star formation efficiency and the fraction of NIRB 
from Pop III stars obtained by imposing that metals ejected by 
Pop III SNe do not exceed the mean observed values in the Lya 
forest. The cut model for the galaxy counts is assumed. 



5 DISCUSSION 

Pop III stars can account for the entire NIRB excess if the 
high redshift star formation efficiency is f+ — 10% — 50%, 
depending on the IMF and on the unknown 'normal' galaxy 
contribution in the L band. Our best results are obtained 
by considering a value z end w 8.8 for the redshift at which 
the formation of Pop III stars ends, almost independently 
of the IMF and of the 'normal' galaxy contribution model. 
A hard upper limit z en£ j < 9 is set by the J band data. The 
above efficiency values are not at odd with expectations for 
objects forming in the Dark Ages (Madau, Ferrara & Rees 
2001), although they seem to be higher than those found at 
lower redshifts. 

In addition, the values of /* found should be considered 
as upper limits since the unaccounted NIRB level could be 
lower due to a different zodiacal light subtraction. To quan- 
tify this uncertainty we have performed fits using the Wright 
zodiacal light model. In this case the DIRBE data alone are 
considered, as NIRS data are only available for the Kelsall 
model. We find that required star formation efficiencies de- 
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crease to 0.15 - 0.41 (z end = 8-10) with a Salpeter IMF, 
0.06 - 0.16 with a heavy IMF, and 0.04 - 0.1 with a heavy 
IMF; at z en£ j = 8.8 the resulting star formation efficiencies 
are 0.23 (Salpeter IMF), 0.09 (heavy IMF), and 0.06 (very 
heavy IMF). However, the errors are doubled with respect 
to the case of Kelsall model for the zodiacal light emission. 
Also, the x 2 analysis does not allow to evaluate a unique 
redshift at which Pop III formation ends. This is due to the 
fact that the break around 1 (im is not as sharp as in the 
previous case, as already mentioned. 

As a consistency test, we have calculated the associated 
IGM metal enrichment. If the whole unaccounted NIRB is 
due to Pop III stars, then we predict that the low density 
('true') IGM would be enriched at the observed mean value 
already by z — 15 — 20, and that the Lya forest data imply 
that either a Salpeter or a heavy IMF with characteristic 
mass M c = 15M are acceptable, provided pair-instability 
SNe are the only source of heavy elements. This conclusion 
can be modified by the two following occurrences: (i) only 
stars with masses in excess of 260 Mq are formed, that lock 
their nucleosynthetic products into the VMBH remnant, ef- 
fectively breaking the light-metal production link (Schneider 
et al. 2002) or (ii) metals are very inhomogeneously mixed in 
the IGM, with filling factors much smaller than unity (Scan- 
napieco, Ferrara & Madau 2002) at these high redshifts. In 
both cases, the limits obtained from metallicity arguments 
are much weaker if not purely indicative. 

To exemplify case (i), let us consider an IMF of the form 
4>{M) = 5(1000 M Q ). By imposing that the whole NIRB is 
due to Pop III stars (/* = 0.04±0.01 and z end = 8.83±0.01) 
and assuming the cut model for the galaxy counts, we obtain 
a density of VMBH Q V mbh(z = 8.8) ~ 7x 10~ 4 . This value 
is two orders of magnitude greater than the measured mass 
density in super massive black holes (SMBHs) found from 
the demography of nuclei of nearby galaxies (Magorrian et 
al. 1998; Gebhardt et al. 2001), Qsmbh = 10~ 4 M7 6 = 2.66x 
10~ 6 (Merritt & Ferrarese 2001), but it is unlikely that all 
VMBHs merge into SMBHs (Schneider et al. 2002; Volonteri, 
Haardt & Madau 2002). Hence it is not clear if this scenario 
can represent a viable solution. 

In the second case (ii), the measured metallicity might 
not be representative of the actual cosmic metal density and 
might result in a large overestimate of the amount of IGM 
heavy elements. In addition, note that above estimates as- 
sume that all the produced heavy elements can escape into 
the IGM. If the metal escape fraction is equal or less than 
~ 10%, the whole unaccounted NIRB could easily be due 
to Pop III stars with no conflict with the observed IGM 
metallicity limits. 

As a final remark, we point out that additional sources 
might contribute to the NIRB, as accretion onto black holes 
(Carr 1994) or decay of massive primordial particles (Bond, 
Carr & Hogan 1986). 
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